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CN ; ABSTRACT 

in 

We estimate the dynamical masses of several young 10 Myr) massive star clusters in two interacting galaxies, NGC 4038/4039 ("The 
QQ Antennae") and NGC 1487, under the assumption of virial equilibrium. These are compared with photometric mass estimates from /f-band 
photometry and assuming a standard Kroupa IMF. The clusters were selected to have near-infrared colors dominated by red supergiants, 
and hence to be old enough to have survived the earliest phases of cluster evolution when the interstellar medium is rapidly swept out from 
. ^ the cluster, supported by the fact that there is no obvious Ho- emission associated with the clusters. All but one of the Antennae clusters 
have dynamical and photometric mass estimates which are within a factor » 2 of one another, implying both that standard IMFs provide 
5_j a good approximation to the IMF of these clusters, and that there is no significant extra-virial motion, as would be expected if they were 
, rapidly dispersing. These results suggest that almost all of the Antennae clusters in our sample have survived the gas removal phase as bound 
or marginally bound objects. Two of the three NGC 1487 clusters studied here have M^,,,, estimates which are significantly larger than the 
photometric mass estimates. At least one of these two clusters, and one in the Antennae, may be actively in the process of dissolving. The 
process of dissolution contributes a component of non-virial motion to the integrated velocity measurements, resulting in an estimated M^,,,, 
which is too high relative to the amount of measured stellar light. The dissolution candidates in both galaxies are amongst the clusters with the 
lowest pressures/densities measured in our sample. 

Key words, star clusters - dynamical masses - NGC 4038/4039 ~ NGC 1487 - IMF 



1. Introduction 

Despite being the targets of intensive studies over the last fif- 
teen years, young extragalactic star clusters, which are found 
in large numbers in interacting galaxies (e.g. 'Holtzman et al.' 
1992^Whitmoreetal. 1993, 1999, Zepf and Ashman 1999, 
Mengeret'al Jl2005i iBastian et al.H2006i iTrancho'et all I 
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65.N-0577, 66.B-0419, 67.B-0504 and 68.B-0530 



as well as in other e nvironments like normal spirals 
(ILarsen & Richtlerl |2004 iLarsen et alj |2004|) . seem to have 



raised more questions than they have answered. 

One of the most obvious, but arguably most interesting, 
questions is: how many of the young star clusters (YSCs) sur- 
vive to old age (i.e. become globular clusters), and what hap- 
pens to the others? Most likely many clusters disperse, con- 
tributing to the general field star population; although, it re- 
mains uncertain what fraction of the general field population 
originated in stellar clusters. 
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Several studies (e.g. lLarsen et al .12004 iBastian et al ]l2006h 
have shown that the properties of (at least some) young clus- 
ters are consistent with them being the progenitors of what we 
see as globular clusters today. Is it possible to identify from a 
population of extragalactic young star clusters those which will 
survive for a Hubble time? Or, as a different way to phrase the 
same problem: How many star clusters, and with which prop- 
erties, formed the host population of the globular clusters in 
today's galaxies? 

In environments as different as mergers like NGC 
4038/4039 and the Milky Way, it seems that, at least up to 
around 100 Myr, 50-90% of the star clusters are destroyed 
within each decade of time. This effect has been named "in- 
fant mortality". 

The current hypothesis ("Hills' 198^, 'Lada et al.' 11984' 
Boilv & Kroupa 2003ab, Fall, Chandar & Whitmore 2005, 
Goodwin and Bastianll2006l I Whitmorel2007h is that the gas re- 
moval caused by stellar winds and supernovae unbinds some 
of the clusters, and that this process is dominant only out to 
roughly 30 Myr. Later, the much slower and less destructive 
two-body relaxation takes over, which has, together with other 
effects like the impact of the galactic gravitational field etc., the 
potential to dissolve another fraction of the initial survivors. 

More observational data are necessary to get a clearer 
idea of the dynamical processes at work during cluster for- 
mation and destruction. The cluster populations analysed so 
far with respect to their ages have not been corrected for 
the (unknown) cluster formation history. However, all studies 
which analyse statistically significant numbers of high-mas s 
clusters (NGC 4038/4039. Pall. Chandar & Whitmorel (l2005b. 
Mengel et all (l2005b . M51, Bastian et al. (! 2005b ) are interact- 
ing systems where the star/cluster formation history is neither 
constant, nor a delta burst, but rather some more complex, un- 
known funtion of time. This certainly affects the age distribu- 
tion of clusters and hence the destruction rate derived from it. 

A different approach targets individual star clusters for in- 
tense studies of their physical parameters, with the goal of us- 
ing these paramters to decide whether a star cluster is doomed 
or a candidate for a future GC. 

2002h 



In our original study (M engel et al 



(as ha d, for example. iHo & Filippenkd ( 1996allbl) . Sternberg 



(119981) ) that clusters are in virial equilibrium, since at ages of 
around 8 Myr, they have survived for many crossing times. 
However, in the view of the high cluster destruction rate de- 
rived from recent studies, this assumption may not be univer- 
sally applicable. 

Our current study expands the number of analysed individ- 
ual clusters. With a larger sample, we hope to be able to find 
a diagnostic to determine the dynamical state of an extragalac- 
tic, and hence only barely resolved, star cluster. Or at least to 
see which potential techniques are unfeasible or do not lead to 
useful results. 

The galaxies we targeted are NGC 4038/4039 and NGC 
1487. While the first, also called "The Antennae", is one of the 
best studied nearby mergers, NGC 1487 is less well known. It 
is described as a peculiar gal axy, which shows t wo faint tails as 
a tracer of earlier interaction. iLee & Led (l2005b conclude from 
their two-colour analysis of the cluster system that the merg- 



ing process could have taken place 500 Myr ago. Most of the 
star clusters are found in three or four "condensations", and the 
brightest clusters, like those targeted for our study, are much 
bluer than the lar ger population of fainter clusters. In total, 
iLee & Led (l2005b found more than 500 cluster candidates in 
HST/WFPC2 data. The galaxy is at approximately half the dis- 
tance to the Antennae, but substantially fainter: Its total magni- 
tude is comparable to the LMC. 

The stellar velocity dispersion in the clusters is typically 
around 15 km s ' and therefore requires medium- to high spec- 
tral resolution of these faint targets, which is only achievable 
with 10m class telescopes. Apart from near-infrared imaging 
for the cluster photometry, we need an estimate of the cluster 
size, which for objects at distances between 10 and 20 Mpc and 
sizes of 2-4 pc requires very high spatial resolution. 

2. Observations & Data Reduction 

In this section we present the ground based imaging and 
spectroscopic data obtained for the K-band bright clusters 
with strong CO absorption (as a consequence, the ages span 
a narrow range around 8.7 Myr) selected for spectroscopy. 
Table [T] lists integration times for both imaging and spec- 
troscopy. Supporting archival images taken with the Hubble 
Space Telescope (HST) are also described. 

2.1. SOFI and ISAAC imaging data 

ISAAC/VLT imaging of NGC 4038/4039 was performed in 
ON/OFF mode during the nights 15.04.2001 (Ks-band) and 
16.04.2001 (CO-bandhead filter). The target fit completely 
onto the detector (0'/1484/pixel, total field size 2^5 x 2.'5). 
Seeing was excellent during both of these photometric nights 
(the FWHM of the PSF from coadded frames is <a.'4). 
SOFI/NTT imaging of NGC 1487 covered J, H, Ks broad- and 
Bry, C02.32jum, PayS, and continuum NB2.28//m, NB1.215//m 
and NB2.195//m narrow bands. Here the field size was roughly 
twice that of ISAAC (4 .'9 x 4.'9, with a pixel size of 0'.'292). 
With the target spanning only a bit more than 2', an efficient 
on-chip offset pattern could be used. 

Reduction of the ISAAC and SOFI broad- and narrow band 
data was performed using the IRAF packag^l]. It included dark 
and sky subtraction (either using the median of several neigh- 
bouring sky images or, where this led to residuals, doing pair- 
wise subtraction), and flat fielding by a normalized median of 
all sky frames. All of the ON frames were slightly offset with 
respect to each other, in order to minimize the effect of pixel 
defects. Therefore, they had to be shifted to a common location 
before using the imcombine task (setting the minmax rejection 
algorithm to reject the highest and the lowest pixel) to combine 
the single frames. 

Photometric standards GSPC S279-F, Ks-magnitude 12.03 
(ISAAC) and S301-D, Ks-magnitude 11.79 (SOFI) from 



' IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Fig.l. False-colour image of NGC 4038/4039 with HST/ACS F814W covering both the blue and the green channel, and 
VLT/ISAAC Ks in the red channel. Clusters which are presented in this publication are marked, using the naming convention as 
in W99 (those which had been Usted there), or according to slit number. 



Persson et al.l (119981) were used for flux calibration of all the 
broad-band data. The resulting zeropoints in Ks (the only ones 
we will give here, because Ks photometry will be used in the 
cluster analysis) were 24.28 mag and 22.27 mag, respectively). 

The target clusters were selected to have a high CO(3-l) 
band-head absorption equivalent width (which is covered by 
the ISAAC and SOFI NB2.34A/m filter), which revealed clus- 
ters at ages 10 Myr, where the near-infrared emission is dom- 
inated by red supergiants. Clusters at ages which are dominated 
by very hot young stars do not show photospheric absorption 
features, or their absorption lines are rotationally broadened, 
making measurement of their stellar velocity dispersion very 
difficult. Locations of the selected clusters are shown in Figs.[T] 
and 111 They must be detected in at least the I-band with HST, 
in order to measure their size (see 93. 2b . 

2.2. Spectroscopic data 

Spectroscopy was performed with ISAAC at VLT-ANTU in 
04/2000, 04/2001 and 12/2001. ISAAC was configured to have 
an C.'S wide slit, and a central wavelength of 2.31/im with a 
total wavelength coverage from 2.25 to 2.37yum. This was suf- 
ficient to include the '^CO (2-0) and (3-1) absorption bands 
at a spectral resolution A/AA~9QQQ. Observations of late-type 



supergiant stars were taken so that they could be used as 
templates for the determination of the velocity dispersion^ 
Observations were performed by nodding along the slit and 
dithering the source position from one exposure to the next. 
B5V atmospheric calibrator stars were observed several times 
during the night. 

The reduction of ISAAC spectroscopy data also made use 
of the IRAF data reduction package. It included dark subtrac- 
tion and flat-fielding (using a normalized flat-field created from 
internal flat observations) on the two-dimensional array. Sky 
subtraction was performed by typically using the median of 3- 
4 frames taken at the other nodding position, and pairwise if the 
first strategy left strong residuals. This was followed by a rejec- 
tion of cosmic ray hits and bad pixels. The spectra were then 
corrected for tilt and sUt curvature by tracing the peak of the 
stellar spatial profile along the dispersion direction and fitting a 
polynomial to the function of displacement versus wavelength. 
Our suite of velocity template stars was used for this purpose. 

Wavelength calibration was performed in a similar manner 
Here, we used a combination of arc discharge lamp observa- 
tions and night sky Unes for the identification of wavelengthts. 



- This research has made use of the SIMBAD database, operated at 
CDS, Strasbourg, France 
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Fig. 2. Left: NGC 1487 false colour image, with HSTAVFPC2 F435W in the blue,HST/WFPC2 F814W in the green, and 
NTT/SOFI Ks in the red channel. HST images were smoothed to the FWHM of the SOFI images. Right: False-colour image 
with the SOFI Ks image in the red, and HST ACS/HRC F814W (not smoothed, therefore substantially better spatial resolution) 
in green/blue. The clusters for which we present spectra are labelled in the image on the right: Clusters 1-3 were observed 
with ISAAC, cluster 1 additionally with UVES. During the largest fraction of the observations, the slit was oriented such that it 
covered all three clusters. Cluster 3 was covered during another integration which accounts for roughly 20% of its exposure time 
(see TabQJ. 



Table 1. Integration times of all images and spectra used for the main analysis of the clusters. We did not list the many different 
narrow band filters which were used for age dating (Bry, C02.32, Pa^S, Ho-). 



Cluster 


Integration times VLT/NTT 


Integration time HST/ACS 




ISAAC spectroscopy 


ISAAC/SOFI Ks imaging 


F814W 


F550M 


F435W F555W 


[W99]2 


2400s 


360s 


3360s 


2544s 


1530s 


[W99]15 


16800s 


360s 


3360s 


2544s 




Sl.l 


28800s 


360s 


3360s 


2544s 




Sl_2 


28800s 


360s 


3360s 


2544s 




S1.3 


28800s 


360s 


3360s 


2544s 




Sl_4 


28800s 


360s 


3360s 


2544s 




2000.1 


16800s 


360s 


3360s 


2544s 




S1.5 


9000s 


360s 


3360s 


2544s 




S2_l 


9000s 


360s 


3360s 


2544s 




S2_2 


9000s 


360s 


3360s 


2544s 




S2_3 


9000s 


360s 


3360s 


2544s 




NGC 1487-1 


14700s 


100s 


640s 




1540s 


NGC 1487-2 


14700s 


100s 


640s 




1540s 


NGC 1487-3 


18300s 


100s 


640s 




1540s 



Single object spectra were combined by shifting-and- 
adding, including a rejection of highest and lowest pixels. 

The object spectra were then extracted from user defined 
apertures (usually the limits of the apertures were placed at the 
point where the counts in the combined spectrum had dropped 
to 1/10 of the peak value). A linear fit to the background (below 



a;7% of the peak intensity for all clusters) on both sides of the 
object spectrum was subtracted. 

For cluster NGC 1487-3, care was required in both the slit 
positioning and the spectrum extraction due to the nearby faint 
companion cluster For one slit orientation, covering the largest 
fraction of the exposure time, clusters 2 and 3 were lined up 
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in the slit, which offset it just slightly from the brightest clus- 
ter 1 (see Figs. |2] and [8]l. For 3600s, we integrated in an al- 
most perpendicular slit orientation which included only cluster 
3. Despite trying to avoid flux from the fainter companion from 
entering the slit or the extraction aperture, and excellent seeing 
during most of the exposure time, it is likely that some un- 
known contribution from the companion is present in our final 
spectrum for cluster NGC 1487-3 (see also Sect. l3.3l l. 

An atmospheric calibrator (B5V) was observed and re- 
duced in the same way as the target and used to divide out the 
atmospheric absorption features from the spectra. 

This article concentrates on ISAAC spectroscopy, and we 
are using only the UVES velocity dispersions previously pub- 
lished dMengel et all I2OO2I l2003i) for clusters where both, 
ISAAC and UVES spectroscopy was perfomed: [W99]-2 and 
NGC 1487-1. Data reduction and analysis of the UVES spectra 
is described in those pubUcations. 

2.3. Hubble Space Telescope Imaging data 

We use imaging observations taken with HS T to estimate the 
size of each cluster. Observations for the size determination 
of the NGC 1487 clusters were obtained as part of Hubble 
Space Telescope Cycle 11 observations (Proposal-ID 9473, 
PI: Vacca). We ai-e using the F814W and F435W Advanced 
Camera for Surveys / High Resolution Channel (ACS/HRC) 
images. For the Antennae clusters we use HST/ACS-WFC im- 
ages obtained for Proposal-ID 10188 (PI: Whitmore) in the 
F550M and F814W filters. In addition, we took advantage of 
higher resolution ACS/HRC images in the F555W filter of a 
supernova in the Antennae (Proposal ID 10187, PI: Smartt), 
which happened to include one of our target clusters ([W99]- 
2). 

For all HS T data, we use the pipeline reduced images. Total 
integration times are listed in Tab[T] The total field coverage for 
the Antennae WFC F814W and F550M images was roughly 
3.'5x3.'5, at a pixel size of O'.'05/pix. The HRC images taken for 
NGC 1487 and the Antennae supernova had a total field size of 
around 31"x31", at a pixel size of 0'/027/pix. 

For p hotometry, we use the p hoto metric zeropoint deter - 
mined by iDe Marchi et alj ( |2004|) and ISirianni et all (l2005b . 
The reduced HS T images were combined with our Ks-band 
images in order to create the two-colour-images shown in 
Figures [T] and |2] 



3. Analysis 

3.1. Velocity Dispersion Measurements 

For each cluster spectrum, we estimated the width <t of the 
broadening of the CO absorption features (assumed to be a 
Gaussian) which best fit the cluster spectrum in the following 
way. An appropriate stellar template spectrum (described be- 
low) was broadened by Gaussian functions of variable cr, rang- 
ing from to a few 100 km/s, and shifted in wavelength by 
radial velocities between 1400 and 1800 km/s. We set the start- 
ing point at 15 km/s and a radial velocity of 1600 km/s. The 
resulting set of broadened templates were then compared with 



X SNR . 50, 2-29-2-325 Mm (CO bandhead 
X SNR = 20, 2-29 -2-325 iim (CO bandhead) 
X SNR = 20, 2-25-2-285 fim (Ca, le, Mg etc) 




Fig. 3. Results of the fits of artificially broadened spectra, 
which had noise added corresponding to signal-to-noise-ratios 
(SNRs) of 50 and 20, respectively (where SNR*20 is the SNR 
obtained for most of our clusters). At our typical S/N, strong 
deviations are only expected at less than half of the resolution 
of our data (or <6-7 km s"'). 



the cluster spectrum. The best fit was determined by evaluating 
X' and then searching for the minimum of the function ;^'^(v,-,cr) 
using a simplex downhill algorithm. Radial velocities are given 
in Table 111 velocity dispersions in Table |4] 

It is important that the template spectrum be a good overall 
match in terms of stellar features to the cluster spectrum. For 
a star cluster that formed ~10 Myrs ago, late K through early 
M supergiant stars are expected to provide the largest contribu- 
tion to the 2.3/jm flux . However population synthesis models 
(iLeitherer et al. 1 11999* ) also show that hot main sequence stars 
will make a non-negligible contribution to the flux at this wave- 
length. Since hot O and B-type stars have an essentially feature- 
less spectrum in the near-infrared, they only represent a dilut- 
ing continuum, which decreases the equivalent width of the CO 
band-heads. This has the eff'ect of shifting the apparent domi- 
nant stellar type towards higher effective temperatures. Starting 
out with a template spectrum with weak CO features leads to 
very low velocity dispersions, while the the opposite is true if 
an M5I star (which has strong band-heads) is used, with sub- 
stantial differences in the results (a few km s"' to up to about 
30kms-'). 

We believe that no significant bias has been introduced 
through the selection of a stellar template and/or the wave- 
length range which was considered, because we have a large 
suite of templates which allowed us to find a good match for 
each cluster. For a good match the best fitting velocity disper- 
sion was essentially independent of the selected wavelength 
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Fig. 4. ISAAC spectrum of NGC 1487-1 with a good fit (using an MSI template spectrum) with a velocity dispersion of cr = 13.7 
km/s as the bottom spectrum, and a bad fit (1 1 km/s) at the top. The residuals are for the good fit. 



range, see iMengel et al.l (l2002h for details. While we cannot 
rule out slight mismatches between template and cluster spec- 
tra, particularly for the low-SNR clusters S2_3 and Sl_5, these 
should only have a slight effect on our velocity dispersion mea- 
surements, and hence Mdy„ estimates. 

The determination of the velocity dispersion from the op- 
tical echelle data used the same procedure as for the ISAAC 
spectra, relying on the Calcium Triplet around 8500A, but also 
using the Mgl absorption feature at 8800A and other weaker 
metal absorption lines between 8400 and 9000A. 

Clusters [W99]-2 and NGC 1487-1 were observed with 
both instruments, ISAAC and UVES, in order to check the 
consistency of the final velocity dispersion measurements. The 
agreement was excellent for the Antennae cluster {cr - 13.9 + 
2.0 km/s from ISAAC and cr = 14.3 ± 0.5 km/s from UVES) 
and reasonable for NGC 1487-1 (cr = 12.0 + 2.0 km/s and 
cr = 15.4 + 2.0 km/s), revealing no strong systematic instru- 
mental eff'ect. 

Our velocity dispersion measurements are summarized in 
Table |4] We find a range from about 7 km s ' to more than 
20 km s For the optical echelle data the instrumental resolu- 
tion is only cr instrument-^ -"^ km s"', ensuring that the cluster line 
profiles are always well-resolved. The situation for the near-IR 
data is somewhat less straightforward, because the instrumental 



-- 14.2 km s , and some of our measure- 



resolution is cr,„ 5,,.,,, 

ments lie below this value. However, as shown in iMengel et al 



(2002) and Fig. [3] results are reliable and reproducible down to 
approximately half this value. 

Fig. 13 supports the assumption that the origin of this ef- 
fect is the large amount of signal contained in the CO band- 
heads. We used our stellar spectra to create artificially broad- 
ened, noisy spectra, and re-determined the velocity dispersion 
100 times for each sampled input velocity dispersion. We used 
two different wavelength ranges of equal length, one including 
the CO bandhead, and one which only includes only several 
metal absorption lines. The fit using the CO region has a stan- 
dard deviation which is ^ five times smaller than that from the 
other region. Mean fit values from the two different wavelength 
regions start separating around cr-\6 km/s, becoming quite 
substantial below 10 km/s. The region around the CO band- 
head allows re-determination of the velocity dispersion down 
to around 6 km/s. Therefore we believe that our measured ve- 
locity dispersions, within their reported uncertainties, are valid, 
because they all lie at or above the limit for a reliable CO fit. 

Figures|4]and|5]show some of the fits to the observed cluster 
spectra. While Fig. H] shows a good fit in the middle (with the 
corresponding residuals at the bottom), and a bad fit at the top. 
Fig. |5] shows the best fits for all our spectra. In general, we 
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2.30 



2.32 



Rest wavelength [jxm] 

Fig. 5. All ISAAC spectra, overplotted with their best fits, with labels. The region shown here generally provided the best fit. 
Table 2. Radial velocities and J2000 coordinates of the Antennae clusters. 



Cluster 


Yr (observed) 


Heliocentric correction 


V, (corrected) 


RA 


Dec 




km/s 


km/s 


km/s 


hms 




[W99]2 


1628.5 ±0.5 


-9.8 


1618.7 


12:01:50.46 


-18:52:13.89 


[W99]15 


1621.3 ±0.3 


-9.0 


1612.3 


12:01:55.43 


-18:52:19.45 


Sl.l 


1604.0 ±0.2 


-8.1 


1595.9 


12:01:55.20 


-18:52:18.35 


SO 


1602.0 ±0.2 


-8.1 


1593.9 


12:01:55.26 


-18:52:20.04 


S1.3 


1594.3 ±0.3 


-8.1 


1586.2 


12:01:55.49 


-18:52:25.49 


S1.4 


1608.0 ±1.2 


-8.1 


1599.9 


12:01:55.05 


-18:52:14.67 


2000.1 


1618.6 ±0.2 


-9.0 


1609.6 


12:01:55.29 


-18:52:17.18 


S2.1 


1580.7 ±0.2 


-8.6 


1572.1 


12:01:55.54 


-18:52:20.73 


S2J1 


1644.5 ±0.4 


-8.6 


1635.9 


12:01:51.90 


-18:52:27.94 


S2.3 


1658.2 ±0.2 


-8.6 


1649.6 


12:01:53.13 


-18:52:25.63 


S1.5 


1645.5 ±2.0 


-8.1 


1637.4 


12:01:54.31 


-18:51:56.86 



obtained good and stable fits for the clusters, even - surprisingly 
- for the low SNR spectra for SI. 5 and S2.3. 

3.2. Cluster Sizes 

Sizes from ACS images were determined using the routine 
ishape, implemented in the data redu ction package baolab, de- 
veloped by S. Larsen dLarsen Il999h . The clusters all appear 
slightly-to-well resolved in the HST images, ishape convolves 
a user-provided PSF with an analytic cluster profile, and de- 



termines the minimum for a range of sizes using a sim- 
plex downhill algorithm, ishape rejects pixels which deviate 
strongly from the median value of pixels at the same radius 
outside a "clean" radius, which we set to 3 pixels. The fit is 
performed out to a radius which we set to 10 pixels for most 
clusters, 12 pixels for [W99]-2 and N1487-1 (which were the 
brightest clusters of each target, therefore the SNR was suffi- 
cient out to a larger radius). We have taken into account the 
ellipticity which was assumed for the best fit model, and de- 
termined the projected half-light radius which would corre- 
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spond to a circular model by using the average of r„„„ and 
r^ax- The general validity of this approach remains to be ver- 
ified, but is currently justified by geometric considerations and 
some numerical integrations (Larsen 2003). We converted the 
output FWHM from ishape to a half-light radius by applying 
the appropriate concentration-dependent conversion factor, as 
described in Larsen (2001). We convert these effective radii 
from arcseconds to parsec by assuming distances of 19.3 Mpc 
(Antennae) and 9.3 Mpc (NGC 1487). 

For the NGC 1487 ACS/HRC data, we created our PSF 
from archival data of a moderately bright star (same filter and 
camera as the science data) which was obtained for a com- 
pletely different purpose (Proposal-ID 10198, PI: Wozniak). 
The Antennae ACS/WFC data contained several foreground 
stars which we used as PSF. Despite careful shifting-and- 
adding, a PSF created from more than two or three stars 
was always slightly broader than the original PSFs. Therefore 
our PSF reference was created from only three stars. For the 
Antennae ACS/HRC data, it was a lucky coincidence that the 
supernova was located in the direct vicinity of cluster [W99]-2, 
thereby providing a suitable PSF reference. 

Our PSFs had the following characteristics: For the 
Antennae images, the FWHM were 2.0 pix (O'/IO) for F814W, 
1.89 pix (a.'095) for F550M, and 2.48 pix (QI!Q66) in the 
F555W filter, where the first two were using the WFC, and the 
third the HRC. Only the F555W PSF showed a slight deviation 
from a Gaussian profile, there was a slight ( 5% of the peak) in- 
crease in the count level ^3.5 pix north of the peak. FWHM of 
the PSFs used for the NGC 1487 images were 2.58 pix (a.'07) 
in the F814W filter, and 1.97 pix (0'/053) in the F435W fil- 
ter. Apart from the strong Airy ring in the F814W PSF, they 
showed no peculiarities. 



Best fits were typically obtained for King (l King|[l962 ) pro- 
files, with concentration parameters between c = 5 and 300. For 
most clusters, we determined the best-fitting half-light radius 
by running ishape for distinct values of c (5, 15, 30, 100, 300), 
and additionally Moffatl5, Moffat25 and Gaussian profiles. We 
did not generally do a two-parameter fit (optimizing c and r/,p 
at the same time), because we believe that this increases the 
risk of getting trapped in a local minimum. But for two clusters 
([W99]-2 and S2_l) we did implement a two-parameter fit, and 
obtained satisfactory results (the concentration for the best fit 
lay between the two best fixed-c fits). 

The projected half-light radii which resulted from the op- 
timization in the two different filters are listed in Table [3] to- 
gether with concentrations. The agreement between the filters 
was exceUent for all the NGC 1487 clusters, [W99]-15 and 
Sl_2, reasonable for [W99]-2 and S2_2, and not very good for 
three clusters, Sl_5, S2_l and S2_3, for reasons that are not 
clear. From visual inspection, the fit in F550M looked much 
better for S2_l, which is why we use this fit as the final value. 
But for the two other clusters, all fits look quite reasonable, and 
we used averages (with rather large uncertainties) as final val- 
ues. Cluster sizes range from x \ - 8pc, with a median size of 
2.9 pc. These values are fairly typical for young s tar clusters 
(e.g. lLarsenl l2004. iLee. Chandar & Whitmorell2005l) . 



3.3. Clusters with multiple components 

For essentially all of our clusters, faint additional point sources 
(or slightly resolved objects) can be detected within the slit 
width of 0'.'3 (see Fig. fTTT l. While in most cases, these objects 
only contribute a few percent of the flux of the primary, we 
have indentified several objects where the impact of the nearby 
cluster is possible but difficult to quantify, or is obvious and 
perhaps severe. Two of the clusters were discarded from fur- 
ther analysis for this reason: Sl_3 and Sl_4 which had been 
only secondary targets in slit position S 1 . The multiple compo- 
nents which show up in the ACS image (see Fig. fTTT l are also 
obvious, even though not as well resolved, in the ISAAC K- 
band image as extended and elongated PSFs. 

This is not the case for the other two Antennae clusters with 
multiple components in the F814W image: For both, Sl_l and 
2000_1, the PSF looks symmetric, and there is no obvious in- 
dication of multiplicity. This could mean that the companion 
clusters are bluer than the main component and therefore their 
contribution in K-band, where we estimate both the velocity 
dispersion and the photometry of the cluster, may be negligible. 
Nevertheless, we have marked these clusters in all our analysis 
plots, even though their properties do not turn out to be unusual 
in any obvious way. 

N 1487-3 is a candidate which is expected to have suffered 
some impact on its measured velocity dispersion: While we 
think that photometry took reasonable account of the multi- 
plicity, and only uses the flux from the primary component to 
estimate the magnitude and photometric mass of the cluster, 
some flux from the companion cluster leaked into the slit dur- 
ing spectroscopy and cannot be removed. This is the case at 
a low level during the major fraction of the integration time 
where the seeing was exceUent (FWHM 0'.'3-0'/4 in K-band), 
and more pronounced during the ^30% of the integration time 
where the seeing was larger (FWHM 0'.'5-0'.'6 in K-band). The 
distance between the two components is 0'.'5. 

We believe that we see the impact of an erroneously high 
velocity dispersion (caused by the different radial velocity of 
the companion cluster) in the very high ratio of M^/jn/Mp/,. 

3.4. Dynamical Cluster Masses 

The dynamical mass of each cluster is estimated using 

TjFcr'^rhp 



M, 



d\n 



(1) 



where rj - rjic) is a factor that depends on the distribution of 
the stellar density with radius (described below), cr is the stellar 
velocity dispersion, r/,p is the projected half-light radius, and G 
is the gravitational constant. The function rj depends on both 
the cluster concentration, c, and on the mass-to-light ratio as a 
function of radius. We assume that it is justified to split these 
two dependencies into two separate parameters: rj - rjic) and 
F = F(t). Figure |9] shows the variation of rj with concentration 
c. To generate the relationship between the coefficient neces- 
sary to estimate the mass of the clusters, 77, and the ratio of 
the tidal radius to the core radius (the concentration parameter, 
c) for a King model (Fig.|9]l, we solved numerically Poisson's 
equation yielding the density distribution for a variety of Wo 
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Fig. 6. Original image, fit and residuals for two clusters (S2_2 and NGC1487-1) for the best and a bad fit each. Pixels deviating 
substantially from the radial median are excluded from the fit and therefore expected to show up in the residuals, for example the 
two fainter objects in the top left corner of S2_2. The better fits are characterised by a smoother distribution at the central location 
of the cluster. For each cluster, the images are shown with the same greyscale range. 
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Fig. 7. VLT-ISAAC Ks images of Antennae clusters included in our spectroscopy slits. The clusters which are marked with 
boxes are those which have an obvious strong neighbour or which consist of multiple objects. This could have an impact on the 
measured velocity dispersion, and the two most obviously aff'ected clusters (Sl_3 and Sl_4) were discarded from further analysis 
for this reason. 
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Table 3. Projected half-light radii of all clusters which do not have very obvious multiple components. Sizes were measured in 
two different filters (where possible), in one additional filter for [W99]-2, and not at all for 2000_1. The latter was assigned the 
same size and concentration as the two clusters in its vicinity (Sl_l and Sl_2), because the three clusters appear comparable in 
size in the K-band images.'' was measured on an ACS/HRC-F555W image. 
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c(F814W) 
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r„„(F550M) 


c(F550M) 
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pc 
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pc 
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30 


0.92 
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300 


0.87 
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[W99]2'' 








6.8+2 
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0.88 






[W99]15 


1.4+0.2 


300 


0.91 


1.5+0.4 


300 


0.85 


1.4+0.2 


300 


Sl_l 


3.6+0.3 


15-300 


0.70 








3.6+0.3 


150 


SO 


3.6+0.5 


300 


0.69 


3.6+0.2 


15-300 


0.74 


3.6+0.4 


300 


2000.1 














3.6+1.0 


300 


S2_l 


2.0+0.6 


300 


0.68 


3.7+0.5 


152 


0.71 


3.7+0.5 


150 


S2J2 


2.0+0.5 


300 


0.75 


2.9+0.4 


100-300 


0.66 


2.5+0.5 


300 


S2_3 


2.7+0.5 


300 


0.68 


3.9+0.3 


30-300 


0.77 


3.0+1.0 


300 


S1.5 
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15-300 


0.90 


0.3+0.2 
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0.78 


0.9+0.7 


15 










r„„(F435W) 
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NGC 1487-1 


2.7+0.3 


30 


0.83 


3.0+1.0 


30 


0.74 


2.8+0.5 


30 


NGC 1487-2 


1.0+0.3 


100-300 


0.75 


1.4+0.2 


5-15 


0.66 


1.2+0.2 


300 


NGC 1487-3 


2.0+0.3 


5-300 


0.70 


2.2+0.4 


15-300 


0.94 


2.1+0.2 


200 



Table 4. Absolute extinction corrected magnitudes (distance moduli of 31.41 and 30.13 for NGC 4038/4039 and NGC 1487, 
respectively), ages, velocity dispersion cr (average value from ISAAC and UVES for [W99]-2 (14.0+0.8 and 14.3+0.5) and NGC 
1487-1 (16.3+1.5 and 15.4+2.0)). Projected half-light radius r/,p. M,i/> is the virial mass determined from e quation[T| with the Icr- 
uncertain ties given i n colum n AM. M^/, is the mass expected from Starburst99 models (version 5.0, 2005, ( Leitherer et al J[l999 



Vazquez & Leithereiil2005h . instantaneous burst with solar metallicity) for a cluster with the given absolute K-band magnitude 
and age. . 
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^-'-0.5 
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1 


8.7+0.3 
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1 4+" 


103!2f 


2000.1 


-16.8+0.3 


^10 


8.5+0.3 
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1 0+" 
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NGC 1487-1 
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0.6 
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(which is the central potential divided by the velocity disper- 
sion and characterises a King model; see e.g., equation 4-131 
from lBinnev & Tremainelll987l) . Integrating the density profile 
with radius provided the mass coefficient, rj. 

For all clusters in NGC 4038/39 and NGC 1487, 77 varied 
between 5.6 and 9.7. 

The factor F = F(t) describes how 77 varies if the mass-to- 
light ratio varies as a function of radius. 

As described in Fleck et al. I (l2006 h their models indicate 
that the first 10 Myrs of cluster evolution - at least for dense 
clusters like those in our sample - result in a steep increase in 
the factor in Eq.lT] and a more gentle increase after that. This 
effect is caused by mass segregation, which leads to a decrease 



in half-light radius, while the total mass and the half-mass ra- 
dius are largely unchanged - i.e. the mass-to-light ratio varies 
with radius. Even though it is expected to depend on several 
parameters how strongly a cluster segregates (density, IMF, up- 
per mass cutoff, initial radius, number of stars), and only the 
density can be determined a priori, we think that it should be 
expected for our very dense, 10 Myr clusters. 

We applied an average factor o f F = TjilrjQ - 1.3 (derived 
from Fig. 14 in Fleck et al. ( 2006h ) to all our clusters. While 
mass segregation is expected theoretically, it should be noted 
that our barely resolved clusters make it impossible to deter- 
mine observationally if they have undergone mass segregation 
or not, because statistics and crowding will cause the highest- 
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Fig. 8. VLT-ISAAC Ks acquisition image of NGC 1487. It 
shows that the two components of Cluster 3 could be well 
resolved. The slit was positioned on the brighter component. 
Nevertheless, some flux from the fainter component entered the 
slit. 
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log c (c = ryro) 

Fig. 9. Variation of the mass coeflicient, //, as a function of the 
cluster concentration, c (c = n/ro). ro is the King radius, and 
r, is the tidal radius of the King model. Our clusters cover the 
range from log(c) - 0.7 - 2.5. 



mass stars to appear mass segregate d in any strongly centrally 
concentrated cluster ( Ascensoll2008 ). 

The dynamical models do not take into account the con- 
tribution of stellar binary orbital motion to the velocity dis- 
persion; however this contribution is expected to be negli- 
gible, due to the large rnasses of the clusters studied here 
jKouwenhoven & de Griisll2007h . 



3.5. Photometric Ages and Masses 

In order to estimate the mass of each cluster based on our imag- 
ing data, we use our K-band images, since these suffer from sig- 
nificantly less extinction than the V band. We performed aper- 
ture photometry in two different ways: 

For NGC 4038/4039, we used a curve-of-growth technique 
to determine the total magnitudes for a large number of clus- 
ters, resorting to magnitudes from fixed (small) apertures plus 
an aperture correction only in cases where the nearest neigh 



bouri ng cluster was closer than a certain limit (IMengel et al 
2005h. 



Since some of our clusters are located in regions of back- 
grounds which show a gradient (see Fig.fTOb. and only a small 
number of clusters needs to be treated, we additionally used 
a more manual approach (which was the only technique used 
for the NGC 1487 clusters): We chose an aperture size indi- 
vidually for each cluster where the signal from the cluster was 
low enough to be comparable to the noise, and selected a back- 
ground region manually to be representative of the background 
expected at the cluster location. Aperture sizes here ranged 



from 1'.'2 to 1'.'8, roughly 3 to 4.4 times the FWHM of the PSF 
in the Antennae images. For NGC 1487, an aperture size of 3" 
was also roughly four times as big as the FWHM of the images. 

For the Antennae clusters, the two techniques gave identical 
results, except for three clusters: Sl_l is fainter by 0.1 mag us- 
ing the manual technique, Sl_2 is 0.2 mag brighter. Differences 
for those two clusters could be expected, because they are lo- 
cated right on the edge of a variation in background intensity. 
S2_l is also brighter by 0.2 mag using the manual approach, 
and here the reason is not obvious. In all three cases, we used 
the manually determined value, but increased the uncertainty 
estimate in the photometry to account for the differences from 
the two methods. 

Since the diffuse light in NGC 1487 is smooth, and none of 
the clusters have confusing neighbouring clusters, we estimate 
that the photometry for clusters NGC 1487-1 and NGC 1487-2 
is reliable. 

For cluster NGC 1487-3, the K-band photometry is compli- 
cated by the fact that there is a second, fainter cluster so close 
by that the two are only marginally separated in our near-IR 
images. The large uncertainty in the photometry is due to con- 
fusion from this neighboring source. To assign a K-band mag- 
nitude to the brighter cluster, we used an aperture that includes 
both clusters and assumed that the relative brightnesses of the 
two clusters is the same in I-band and K-band. The ratio of the 
peak counts in the acquisition image shown in Fig. |8] supports 
this assumption. 

The absolute K band magnitudes range from -13.4 mag to 
-17.4 mag. The absolute extinction corrected K-band magni- 
tudes are converted to Mpi, by comparison with the absolute 
K-band magnitude predicted by a Starburst99 model for a 10^ 
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Fig. 10. Blowup of the densely populated overlap region in the 
Antennae. In order to enhance the background, we chose a log- 
arithmic scaling. The rectangles indicate regions which we se- 
lected manually as a reasonable assumption for the background 
level at the location of the nearby cluster. 



Mq cluster (instantaneous burst, Kroupa IMF, solar metallic- 
ity) of the same age. Resulting photometric masses are listed in 
Table m 

All clusters are fairly massive, at least compared to young 
star clusters in the Milky Way. For interacting galaxies, these 
luminosities, corresponding to photometric masses between 
8xlO'*M0 and 4.5x10^Mq, are not unusual. Furthermore, 
we selected some of the most massive clusters specifically, 
because, at comparable ages, they are the most luminous 
and therefore most easily accessible to high resolution spec- 
troscopy. 

The fact that the clusters in NGC 1487 are generally fainter 
and less massive than those in the Antennae is expected from 
simple statistical considerations. NGC 1487 has a relatively 
low number of clusters compared to the Antennae - only 1- 
10%. Moreover, because we targeted clusters that are amongst 
the brightest in each galaxy, the likelihood that a system like 
NGC 1487, produces one or more clusters which are compara- 
ble to th e brightest clusters in the Antennae, is comparatively 
low (e.g. lWhitmore et al.ll2000llWhitmorell2007h . 

Extinction turned out to be relevant, at least in K-band, 
only for two clusters, Sl_l (Ay=4.6mag) and 2000_1 (Ay ~ 
10 mag). For all the other clusters. Ay is below 2 mag (with the 
uncertainties AAy < 0.5 mag), which translates to A^- below 
«i 0.2mag, making potential uncertainties of AAa: < 0.05 mag 
small in relation to the overall photometric uncertainty. 

We estimate the age of each cluster by comparing the 
broadband filter measurements, the CO(3-l) bandhead, and 
Bry and Calcium Triplet (CaT, from UVES spectroscopy) 



e quivalent widths vyith the predictions of the Starburst99 
( Leitherer et al.l 19991 Vazquez & Leithere J2005 ) evolutionary 
synthesis models. 

For all clusters, we determined a best fit age between 8 and 
9 Myrs. For clusters with ages a;8.5 Myr, an alternative solu- 
tion exists at 10.5 Myr for all cluster properties which are influ- 
enced by red supergiants (M^,Wco,Wcaj-, V-K, etc.) because 
the quantities are double valued at around this age. Since M;^ is 
the same for both ages, all our conclusions remain unchanged, 
because Mp/,, L^r/M^yn etc. are unaflBCted by this uncertainty in 
age. 



4. Discussion and Implications 

An underlying assumption of any mass estimate based on the 
velocity dispersion is that the cluster is self-gravitating (i.e., 
bound). The age of our clusters, w 10 Myr, is much older 
than their estimated crossing times (tcross ~ rhplcr) of 1-few 
xlO^ yrs, indicating that they have already survived for 20- 
50 crossing times. However, recent results suggest that a large 
fraction of clusters becomes unbound and disperses within 
~ 10 -20 Myr due to the removal of interstellar material; there- 
fore it remains possible that the clusters studied here may not 
be in virial equilibrium. If cluster stars are dispersing, then this 
extra-virial motion will lead to a measured velocity dispersion 
which is higher than would be measured for a bound cluster of 
similar mass. 

One way to assess whether our clusters are gravitation- 
ally bound or show evidence for non-virial motion is to com- 
pare the Lk /M determined from velocity dispersion measure- 
ments with those predicted by population synthesis models. In 
essence, this is a comparison of the dynamical and photomet- 
ric masses. A cluster which has L^/M (based on dynamical 
measurements) which is lower than the photometric estimates 
(from the age of the cluster and its measured stellar light), can 
point to non-virial motion resulting from an expanding, dis- 
solving cluster In Figure [12] we compare our estimated clus- 
ter ages and L/f/M for the Antennae and NGC 1487 clusters, 
with the predictions for an instantaneous burst, solar metallicity 
model from Starburst99 (Leitherer et al. 1999). We show pre- 
dictions for L^/M assuming two different IMFs: a Kroupa IMF 
(solid line) and a Salpeter IMF with 0.1 Mq and 100 Mq lower 
and upper mass cutoffs respectively (dashed line). Within the 
measurement uncertainties, the measured properties agree with 
the model predictions for all but two clusters in NGC 1487 and 
one cluster in NGC 4038/4039. This would stifl be the case if 
we had assumed the distance to the Anten nae which was deter- 
mined from the tip of the red giant branch (ISaviane et al .l2008h 
to be substantially lower than our assumed value, 13.3 rather 
than 19.3 Mpc. Since this lower distance would affect both es- 
timates (lowering Mp/, by a factor 2, and lowering the cluster 
sizes and hence M^,.,, by a factor 1.45), the net effect would 
be that the ratio of M^y„/Mp/, needs to be corrected by a factor 
1.38 (correspondingly decreasing W/M by a factor 0.73). In 
general, this would still lead to a good correspondence between 
the photometric and dynamical estimates, and would leave the 
conclusion the same: 



Mengel et al.: Young star clusters in interacting galaxies - NGC 1487 and NGC 4038/4039 



13 



NGC 4038/4039 

t 

Sl_2 



1 
















t 


* 










Sl_3 




Sl_4 


Sl_5 


S2_l 






m 




i 


S2_2 




S2_3 


2000_l _ A 


1.5" A E 



• 

N14S74 N1487-2 



Fig. 11. HST-ACS/WFC F814W images of the clusters used for our analysis. The clusters which are marked with boxes are 
those which have an obvious strong companion or which consist of multiple objects. This could have an impact on the measured 
velocity dispersion. 
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The good agreement between evolutionary synthesis mod- 
els applied to our clusters, in comparison with dynamical 
masses, suggests that there is no strong variation in the IMF for 
all but three clusters in our sample, and that they have likely 
survived the gas removal phase as bound stellar systems. 

The two NGC 1487 clusters (red stars) and one Antennae 
cluster (S2_l) which are offset below the model predictions 
have dynamical mass estimates which are significantly higher 
than the photometric ones. One possible explanation is that 
these clusters have an IMF which is significantly steeper than 
Kroupa/Salpeter. There is little direct evidence for such an in- 



terpretation, and we believe that it is much more likely that 
these are clusters caught in the act of dissolving. The efficiency 
with which a cluster forms stars will impact the probability that 
it survives the expulsion of its natal gas. For example, clusters 
which form stars at lower efficiencies end up with fewer bound 
stars (i.e. shallower potential wells) relative to the left-over gas 
from formation. Such clusters, as momentum input from the 
massive stars expels the gas, have a lower probability of re- 
maining bound tha n a cluster that formed more stars and had 
less remaining gas. Goodwin and BastianI (|2006) explored the 
connection between star formation efficiency and cluster dis- 
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solution by simulating the N-body dynamics of a cluster after 
the expulsion of gas. Using their results (Bastian, priv. comm.), 
we plot the Li^jM ratios for a Kroupa IMF, solar metallicity 
(Starburst99 models) as the red dotted lines in Figure[T2]for the 
following effective star formation efficiences (eSFEs, defined 
as a measure of how far the cluster is out of virial equilibrium 
after gas expulsion): 60%, 50%, 40%, 30%, 20%, and 10%, 
starting from the top. eSFEs w 40% and higher are predicted to 
result in stable clusters after 20 - 30 Myr, even though many 
clusters, at least the three between the 40% and the 60% lines, 
may lo se a substantial amount of mass dOoodwin and BastianI 
(l2006h'). 



Baumgardt & Kroupal (120071) ran a grid of models with 



larger parameter space, varying star formation efficiency (SFE, 
defined in the normal way as the ratio of stellar mass over 
mass of stars and gas), gas expulsion time and tidal field- 
Even t hough a direct comparison to the lOoodwin and Bastian 
(l2006h results is difficult, some general conclusions are the 
same in both models: All clusters, even the "survivors", 
expand initially, and almost indistinguishably. After 10-20 
Myr, the dissolving clusters continue to expand, while those 
with a sufficiently high SFE re-contract. A more gradual 
gas expulsion than the inst antaneous expulsion assumed by 
Goodwin and Bastian (|200^ makes it easier to remain bound. 

Two clusters in in NGC 1487 and one in the Antennae, in 
Fig. [12] lie in or very c lose to regions where cluster dissolu- 
tion is expected from the lGoodwin and Bastian (12006 ') models. 
We consider two of them candidate dissolving clusters (which 
is particularly interesting because of their high mass - even 
though infant mortality is claimed to be mass independent, 
high-mass clusters are usually intuitively considered more sta- 
ble against dissolution). Variations in the initial conditions, for 
example long er gas expulsion ti mes or tidal fields, as explored 
by Baumgard t & Kroupal (l2007h . would shift the boundary be- 
tween dissolving and surviving clusters in this plot down or up, 
respectively. 

The third cluster is NGC 1487-3, which is likely to have 
suffered effects of cluster multiplicity (see Sect. I3.3l l. 

One of the main results of this work is that most of the 
clusters in our sample appear to have survived, as bound stellar 
systems, the gas removal phase which occurs during the life of 
every cluster It is important to note that, in light of many recent 
works which show that many or most clusters (roughly 50% 
to 90%) probably do not survive the earliest phases of evolu- 
tion, our study targets clusters which are likely to have survived 
this phase. After 10 Myr other mechanisms will continue to 
unbind clusters. If essentially all clusters which reach an age 
of X! 10 Myr in the Antennae are marginially bound or bound 
at this point, this would im ply a very large number of young 
globular clusters. However. iFall. Chandar & Whitmore (l2005() 
show that, at least statistically, star clusters continue to get dis- 
rupted approximately independent of mass, out to an age of 
^100 Myr 

Evaporation of stars resulting from two-body relaxation 
will eventually disrupt a number of lower mass clusters over 
a Hubble time. Clusters with current masses > few x 10^ Mq 
would likely survive this process over a Hubble time, assuming 
the typical evaporation rate of //gy = 1 - 2 x 10"^ Mq yr"'. Such 
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Fig. 12. This is a modifi ed version of a figure from 
Goodwin and Bastian (l2006l) : their Figure 5. It shows L/M for 
a Starburst99 (v.5.0, 2005) SSP with a Ki'oupa IMF, instanta- 
neous burst, solar metallicity (solid, black). Dotted (red) take 
into account the Goodwin & Bastian model results (Bastian, 
priv. comm.) for clusters with effective star formation efficien- 
cies (eSFEs) below 100% (lowest line: 10%, top: 60%). The 
locations of our clusters are indicated (clusters with identical 
ages have been offset by 0.1 Myr in order to be able to see 
which error bar belongs to which cluster). Those which end up 
between the lines according to eSFEs of 20% and 50% are dis- 
cussed as dissolution candidates in the text, and correspond to 
ratios of M,;y„/M/,/, > 2.5. 



a rate is plausible as it reproduces the observed turnover in the 
mass function of globular star clusters in many galaxies (e.g.. 



Fall & Zhang|l200lllWaters et al.ll2006lljordan et alJEooel) . 

Spectroscopic studies like the one presented here require 
enormous amounts of telescope time, and still result only in 
very few spectra. It would be much more efficient if a cluster 
population could be separated into "survivors" and "dissolvers" 
from (high resolution) imaging alone, because this would give a 
better handle on the infant mortality rate. Obviously, constrain- 
ing the infant mortality rate is of immense relevance for the 
whole issue of star formation, because the currently cited clus- 
ter destruction rates range in impact from "cluster formation is 
an interesting, but not very important mode of star formation" 
(for destruction rates of a few tens of percent) to "essentially 
all stars formed in clusters" (for destruction rates of a;90% per 
decade). 

One expectation for expanding, unbound clusters is that as 
they expand, their internal density should decrease. Note, how- 
ever that it is impossible to determine from the size and/or 
concentration of these young clusters alone if they are bound 
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or dissolving: Even though unbound stars, leaving the clus- 
ter with escape velocities of tens of km/s reach distances of 
tens or hundreds of pc in a few Myr, the half-light radius of 
the cluster is not immediately affected severely - even after 
20 initial crossing times, the half-mass radius of a dissolv- 
ing cluster is only 40 % large r than that of a surviving cluster 
dBau mgardt & Kroupa 2007|) . An alternative explanation for 
clusters having low density is simply that they formed in a 
low-density environment, which is expected to lead to smaller 
SFEs. In any case, we might expect our clusters with non-virial 
motions to have low stellar densities as well. In Fig. [13] we 
show the estimated half-mass density p/, (with r/y, and Mpi, as 
input) for the clusters versus the ratio dynamical to photometric 
mass. This figure shows that our two dissolving cluster candi- 
dates (i.e. those with high ratios of MdmlMphot) also appear to 
have low stellar densities. 

We have included in Fig. [13] all the data from the litera- 
ture where the three parameters M^/,,,,, M^/, (assuming a Kroupa 
IMF) and r/,p were provided or could easily be deduced. Most 
of the clusters in other publications are considerably older than 
10 Myr, therefore it is not surprising that they do not show in- 
dication of cluster expansion. But the two literature clusters 
where the dynamical mass exceeds the photometric mass by 
more than a factor two (NGC 6946-1447 and NGC 5236-805) 
confirm the trend shown in our clusters, since they also have 
low densities. 

Clusters with high M^yn/Mp/, likely have low values of p/, 
because they are expanding and thus the high values are due 
to dynamical evolution. Therefore clusters with ages around 10 
Myr which have low densities are excellent candidates for clus- 
ters in the process of dissolving, although clearly some fraction 
of low density clusters appear to be bound at this age as well. 
Low density clusters which are gravitationally bound, such as 
found in the Milky Way ( outer globular clusters , HarrisI ( 19961) ). 
the Magellanic Clouds (Ivan den Bergh 1991), and in nearby 
spirals and lenticular galaxies ( Chandar, Whitmore & Lee 
2004 iLarsen & Brodiell2000l jPenget aL ,20061) . may survive 
longer than their higher density counterparts, since they are ex- 
pected to have lower rates of relaxation-driven stellar evapo 



ration (IMcLaughlin & I2OO7I Ichandar. Fall & McLaughlinI 
2007h . 

A similar picture emer ges if we consider pr essure in- 
stead of density: Following lElmegreen et al. I (l2000h . we used 
P cc GM^i^/rf^^ as an estimate for the pressure in the ambi- 
ent medium during cluster form ation. Then SFE is e xpected to 
also scale with this parameter ( Elmegreen et al J EOOO ). Indeed, 
as shown in Figure[T4] the clusters which may be disolving are 
amongst those with the lowest pressures in both our sample and 
the sample taken from the literature. 

Whereas low pressure/density is not a unique identifier 
for dissolving cluster (since, as shown, there exist also low- 
density/pressure clusters with no sign of expansion), about 50% 
of the clusters in our plots below a density/pressure limit (10^^ 
and 10^, respectively) are dissolution candidates. They consti- 
tute only 20-25% of the whole sample. 

In our future work, we will compare dynamical and photo- 
metric mass estimates of a larger sample, which would provide 
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Fig. 13. The ratio of dynamical mass over photometric 
mass vs. half mass density p/, (taken as the log on the 
X-axis). High values (> 2.5) of Mrf,,,, / Mp/i, which may 
indicate dissolving clusters, seem more common at low 
densities. Symbols are described in the legend, diamonds 
and hexagons represent our data (with the open circles 
marking clusters which may have been influenced by nearby 
clusters), squares are taken from th e literature: M82-F: 
McCr adv. Graham and Vaccal (l2005h . iMcCradv & Graham 



(2007D, M82-A1 



respectively): iSmith et al 



(l2006l) 
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McCradv & Grah^nl (120 07). M82 MGG-9 and MGG-11: 
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320041) 



4214-13, NGC 4449-27 and NGC 
tOOA NGC 5236-502, 5236-805: 



ILarsen & Richtleii (120041) . NGC 6946-1447: 
(l2004l2006l) .NGC 1569-B: lLarsen et al.1 (l2008h . 
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a more robust estimate of the fraction of clusters which appear 
as single entities, but are unbound at an age of xilO Myr. 
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